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BIFURCATION  IN  THE  DUFFING  EQUATION 
WITH  INDEPENDENT  PARAMETERS,  II 


Snopsis : In  a previous  paper,  the  authors  gave  a complete 

description  of  the  number  of  even  harmonic  solutions  of  Duffing's 
equation  without  damping  for  the  parameters  varying  in  a full 
neighborhood  of  the  origin  in  the  parameter  space.  In  this  paper, 
the  analysis  is  extended  to  the  case  of  an  independent  small 
damping  term.  It  is  also  shown  that  all  solutions  of  the  undamped 
equation  are  even  functions  of  time. 

1.  Introduction. 

Consider  the  Duffing  equation  with  damping, 

2 

d u du  3 . , , , , 

— 2 + u = Pxu  + P4  3t  “ P2U  + p3  Cost  (1,1) 

dt 

where  p = Cp^,P2»P3»P4)  is  a real  four  dimensional  vector  varying 
in  a neighborhood  U of  the  origin.  Our  objective  is  to  discuss 
the  2iT-periodic  solutions  of  Equation  (1.1)  for  each  p in  a 
sufficiently  small  neighborhood  U.  In  a previous  paper  [2] , the 
authors  gave  a complete  description  of  the  number  of  solutions 
for  the  undamped  case,  p^  = 0,  under  the  hypothesis  that  the  solu- 
tions were  continuous  in  p^  and  even  in  t.  The  hypothesis  of 
evenness  made  it  possible  to  reduce  the  discussion  to  a single  bi- 
furcation equation.  When  the  damping  term  is  present,  such  an 
hypothesis  is  meaningless  and,  therefore,  two  bifurcation  equations 
must  be  considered.  A complete  analysis  is  possible  after  one  has 
exploited  the  symmetry  properties  in  Equation  (1.1)  to  obtain 
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detailed  qualitative  properties  of  the  bifurcation  equation.  This 
latter  information  is  also  used  to  show  that  all  2ir-periodic  solu- 
tions of  the  undamped  equation  are  even  in  t.  The  analysis  in 
the  present  paper  is  in  the  same  spirit  as  in  [2]  and  uses  ideas 
from  [1] , [3] . 

2 . The  Bifurcation  Equations. 

Consider  the  system: 

. 2 3 

— y + u - p.u  + p_u  = 0.  (2.1) 

dr  X 

A 2TT-periodic  solution  continuous  in  p3,p4  is  a continuous 

function  from  a deleted  neighborhood  V -{(0,0)}  (depending  on 
2 

P1,P2)  of  0 e ^ into  the  space  of  2ir-periodic  functions  which 
associates  to  each  (p^,p2)  e V - {(0,0)}  a 2iT-periodic  solution 
u (P1,P2,P3,P4)  (t)  of  Equation  (1.1).  Furthermore,  the  set 
{u(p1,p2,p3,p4) , (p3,p4)  e V - {(0,0)}},  with  the  uniform  topology, 

is  precompact  and  every  limit  point  of  this  set  as  (p3,p4)  0 

is  a 2TT-periodic  solution  of  (2.1). 

The  idea  for  considering  this  particular  type  of  continuous 
dependence  on  the  parameters  came  for  a paper  of  Hale  and  Taboas  [4] 
in  the  consideration  of  2tt -periodic  solutions  of  another  type  of 
equation.  This  definition  is  more  general  than  the  one  considered 
by  the  authors  in  [2]. 

Since  we  are  only  interested  in  solutions  of  this  form,  it  is 
necessary  to  discuss  some  detailed  properties  of  the  solutions  of 
(2.1)  . 


The  ncessary  information  is  contained  in  the  following  lemma, 
whose  proof  was  given  in  [2] . 

Lemma  2.1.  There  is  a constant  k > 0 and  a neighborhood  U in 
2 

ga  of  Cp1fP2)  = (0,0)  such  that  a nonconstant  2ir-periodic  solu- 
tion u of  (2.1)  exists  for  (p^p^  e U if  and  only  if  pxP2  > 0, 
this  solution  is  unique  except  for  a phase  shift  and  satisfies 


u (t)  | < k|p1/p2 


If  (P1rP2)  e u and  either  p1p2  < 0,  p2  > 0 or  Px  f P2  * 
the  only  2TT-periodic  solution  of  (2.1)  u = 0.  I_f  P-^P2  ? 

p2  < 0,  there  are  the  2Tr-periodic  solution  u = 0, 
u = ±[  (l-P1)/|p2| 11/2.  If  Px  = 0,  p2  = 0,  every  solution  of 
C2.1)  is  27T-periodic. 

For  p2  < 0,  there  are  always  two  27r-periodic  solutions  of 

Equation  (1.1)  which  exist  for  p = (P1»P2»P3»P4)  e ^ in  a 

1/2 

neighborhood  of  zero  and  coincide  with  ± [ (1— p^/ | p2 | ] for 

p3  = p^  = 0.  This  is  proved  in  the  same  manner  as  in  [2]  and  is 
stated  precisely  in  the  following  lemma. 

Lemma  2.2.  There  is  a neighborhood  U in  ^ o£  p = 0 and 
that  for  every  p e U,  p2  < 0,  there  exist  two  27r-periodic  solutions 
of  Equation  (1.1)  which  are  continuous  in  (p3,p4)  and  coincide 
with  ± I (l-PjL) / |P2  | ] 1/2  for  p3  = p4  = 0.  All  other  such  2 it- 
periodic  solutions  of  Equation  (1.1)  for  p3  = p4  = 0 must 
coincide  with  the  solution  u = 0 of  Equation  (2.1)  or  a nonconstant 


I 


I 


periodic  solution  of  Equation  (2.1). 

To  obtain  the  other  27r-periodic  solutions  of  Equation  (1.1) 

-1/2 

continuous  in  (p3,p4),  let  p2  > 0,  u = vp2  in  Equation  (2.1) 

to  obtain  the  equation 


~2  + v = P-jV  + P4  dt  “ V + CT  cost 


(2.2) 


1/2  1/2 
where  a = P2  p^.  From  Lemma  2.1,  we  know  that  p2  u is 

bounded  and,  furthermore,  that  the  only  2-rr-periodic  solutions  of 

Equation  (2.2)  that  need  to  be  considered  are  those  for  which  v 

is  small.  If  p2  < 0,  the  same  remark  is  true  (of  course,  with 

-v^  replaced  by  -t-v^)  provided  the  two  solutions  in  Lemma  2.2 

are  excluded  from  the  discussion. 

We  now  discuss  2iT-periodic  solutions  of  Equation  (2.2)  for 

(v,p^,a,p^)  in  a small  neighborhood  of  the  origin.  The  procedure 

will  be  the  classical  one  used  in  [2]  to  obtain  the  bifurcation 

equations.  Any  2iT-periodic  solution  of  Equation  (2.2)  for  p^  = 

2 

P4  = a = 0 must  be  equal  to  r cos(t-<f>)  + 0 (r  ) for  some  con- 
stants r , <)> . Therefore,  by  letting  t t + <j»,  we  will  obtain 
a solution  of  our  problem  by  considering  the  2n-periodic  solutions 
of  the  equation 


+ v = p.v  + P4  hT  “ + 0 cost(t+<fl)  (2.3) 

dt 

2 

which  for  p1  = p4  = c = 0 are  equal  to  r cost  + 0 (r  ) . 

Let  & = {h:  & : h is  continuous  h(t+2ir)  = h(t)}  and 


*!*»».* gWfci-injtv. 


for  any  h e 9,  let  |h|  = sup|h(t)  |.  Let  P:  & -*■  & be  the 

t 

projection  defined  by 


1 f2lT 

(Ph)  (t)  = — cost  h(s) 

) o 


cos  s ds  + 


1 . f 2l 

— sint 
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h(s)sin  s ds.  (2.4) 


For  any  h e the  equation 


+ v = h 


(2.5) 


has  a solution  in  & if  and  only  if  Ph  = 0.  Furthermore,  there 
is  a continuous  linear  operator  K:  (I-P)i?  -*■  (I-P)i^  such  that 

K(I-P)h  is  the  unique  solution  of 


+ v = (I-P)h 


(2.6) 


which  satisfies  PK(I-P)h  = 0;  that  is,  K(I-P)h  is  simply  the 
2ir-periodic  solution  of  Equation  (2.5)  which  does  not  contain  cost, 
sint  in  its  Fourier  series.  To  this  solution  K(I-P)h,  one  can 
add  an  arbitrary  linear  combination  of  sint  and  cost  to  obtain  the 
general  solution  of  Equation  (2.6).  As  remarked  earlier,  it  is 
only  necessary  for  us  to  add  a term  r cost. 

If  r is  fixed  and  we  define 


P_  = {h  e P:  (Ph)  (t)  = r cost}, 
r 

c t \ . dv  3 

f (v,Pl,p4)  = Pxv  + P4  3^  " v 


then  v is  a solution  of  Equation  (2.3)  in  Pr  if  and  only  if 


v = r cost  + w,  we  (I-P)^ 


+ w = (I-P)f(r  cos  ( • ) + w,plfp4) 


P[f(r  cos(*)  + WfP^p^j)  + a cos(*+<j>)]  = 0 


(2.7a) 


(2.7b) 


(2.7c) 


since  (I-P)  cos  (•+<£)  = 0. 


By  an  application  of  the  Implicit  Function  Theorem,  there  are 
6 > 0,  e > 0,  such  that,  for  |r|  < 6,  | p1 | + | p4 | < e,  there  is 
a unique  solution  w*(r,p^»P4)  of  Equation  (2.7b)  in  (I-P)i^, 
the  function  w*(r,p1,p4)  is  analytic  in  r,p1,p4  and 
w*(0,plfp4)  = 0.  Furthermore,  it  is  very  easy  to  see  that 
w*(r,p^,0)  is  an  even  function  of  t since  v = r cost  + w and 
only  powers  of  v occur  on  the  right  hand  side  of  Equation  (2.7b) 
for  p4  = 0. 

Since  w*(r,p^,p4)  is  uniquely  determined  it  follows  that  there 
is  a solution  v = r cost  + w e P^  of  Equation  (2.3)  which  lies  in 
a sufficiently  small  neighborhood  of  zero  if  and  only  if 
v = r cost  + w*(r,p1,p4)  and  (r,c(>  ,Plfa,p4)  satisfy  the  bifurca- 
tion equations 


[f(r  cos  ( * ) + w*  (r,p^,p4)  , P1»P4)  + cos  ( • + <J> ) ] = 0. 


From  the  definition  of  P in  Equation  (2.4),  these  latter  equations 
are  equivalent  to  the  system  of  equations 


G;LCr,<(),p;L,a,p4)  = a cos  <}>  + — j f(r  cost  + w (r,p1,p4)  (t)  ,plfp4)cost  dt 
G2Cr,<(),p1,afp4)  = cr  sin  <J)  + ^ J f(r  cost  + w*(r,plfp4)  (t)  ,Plfp4)sint  dt 


Since  w*(r,plf0)  is  an  even  function,  it  is  clear  that 
G2  (r , <J> , px , ct , 0 ) = a sin  <J>.  For  w*(r,p1,p4)  = 0,  it  is  easy  to 
evaluate  the  above  integrals.  If  this  computation  is  made  and 
one  uses  the  Taylor  series  to  obtain  the  order  estimate  0(|p1r| 
|p4r|  + | r | 3 ) on  w*(r,plfp4),  the  bifurcation  equations  become 


G1Cr,()),p1,a,P4)  = a cos  <f>  + PjT  “ 4 r + rg1(r,p1,p4)  = 0 
G2Cr,<(),p1,a,p4)  = a sin  $ + P4r  + P4rg2  (r,prp4)  = 0 


where 


gl(r,pl,p4)  = 0 C I ^1 1 2 + IP4 1 2 + I PXP4 1 + r2 1 Px I + r2|P4l  + r4) 


g2Cr,P;L,p4)  = 0(r  + |Pl|  + |p4|). 


These  results  are  summarized  in  the  following  lemma. 


Lemma  2.3.  There  is  a neighborhood  U C W o£  (r,p1,a,p4)  = 
and  a neighborhood  V C of  v = 0 such  that  Equation  (2.3) 


has  a solution  v e Pr  n V for  (r,plfa,p4)  e U and  a given  <j» 
if  and  only  if  (r  ,p^  ,a,p4  , <)» ) satisfy  the  bifurcation  Equations 
(2.8)  where  g. ,g_  satisfy  (2.9). 
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An  immediate  corollary  is  the  following  result  on  the  undamped 
Duffing  equation. 


2 

Corollary  2.1.  There  is  a neighborhood  U C of  (plf  o)  = 0 

and  a neighborhood  V C 9 of^  v = 0 such  that  the  only  2ir-periodic 
& solutions  in  V of  the  undamped  Duffing  equation 


~2  + v = p^v  - v + a cost 


(2.10) 


are  even  functions  of  t if  a ? 0. 


Proof.  For  P4  = 0,  the  second  bifurcation  Equation  (2.8b)  is 
a sin  <j>  = 0.  Therefore,  if  a ^ 0,  we  must  have  <J>  = 0.  If 
<p  = 0,  then  v ( t ) = r cost  + w*  (r,plf  0)  (t)  is  an  even  solution  of 
Equation  (2.10).  Since  all  solutions  of  Equation  (2.10)  can  be 
obtained  from  Equation  (2.3)  and  the  Lemma  2.3,  we  have  proved  the 
desired  result. 


Remark  2.1.  The  conclusion  of  Corollary  2.1  is  also  valid  for  the 
equation 


~ J + V = P-,  v - v3  + of  (t) 
dt 


where  f(t)  is  even  in  t and  f(t)  cost  dt  = tt.  To  prove 

J 0 

this,  one  considers  the  same  equation  with  t replaced  by  t + <t> 
and  observes  that  the  bifurcation  equation  obtained  by  projecting 
onto  sint  has  the  form  a(sin  <j))h(r,4>,p1)  = 0 with  h(0,<)>,0)  = 1. 


i 


-s- 


One  can  now  repeat  the  same  argument  as  in  the  proof  of  Corollary 
2.1  to  complete  the  proof. 


3 . Analysis  of  the  Bifurcation  Equations. 

In  this  section,  we  analyze  the  bifurcation  Equations  (2.8) 
for  (r,4>)  and  determine  the  surfaces  (bifurcation  surfaces)  in 
the  parameter  space  across  which  the  number  of  solutions  (r,<}>) 
changes.  If  the  number  of  solutions  is  to  change  at  some  point 
in  the  parameter  space,  then  there  must  be  a multiple  root  (r,<f>) 
of  Equations  (2.8)  at  these  values  of  the  parameters;  that  is, 

3 (G^,G2)/9  (r,<(>)  = 0.  Our  objective,  therefore  is  to  obtain  the 
values  of  the  parameters  at  which  multiple  solutions  of  Equations 
(2.8)  occur.  Basic  to  the  investigation  is  the  following  a priori 
estimate . 


Lemma  3.1.  There  is  a neighborhood  U C ^ of  the  point 
(r,p^,cr,p^)  = (0,0, 0,0)  and  a constant  c such  that  all  solutions 
(r,<p)  of  Equations  (2.8)  with  (r,p^ , a , p4 ) e U satisfy 


r|  < cC|p1|1/2  + I cr  | 1/3  + | P4  | 1/2 ) 


Proof . If  this  is  not  the  case,  then  there  is  a sequence 
Wn'rn'Pin'°n'Pin>  with  (rn ' pln ' an ' p4n > " (0'°'°'0)  as  n - 


ltn'  n rln  n r4n 
such  that 


- 0, 


"*■  0, 


/ 
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as  n -►  °°.  Dividing  the  Equation  (2.8a)  by  r3  and  letting  n -*■  °°, 
this  implies  that  0 = - , which  is  a contradiction.  This  proves 

the  lemma. 


Lemma  3.1  justifies  the  scaling 


6a 


2/3 


in  Equations  (2.8).  Making  this  change  of  variables  and  dividing 
the  equation  by  a,  one  obtains  the  equivalent  equations  for  a / 0. 


F def  3 3 , r\  t 2/3, 

= pp  - j p + cos  <)>  + p0(a  ) = 0 

G2  dlf  6p  + sin  <{)  + 6p(a2/3)  = 0. 


For  (p,a,6)  = (0,0,0),  these  equations  only  have  the  solution 
(p,4>)  = ( (j) 1//3 , 0)  . Furthermore, 


3 (G^ , G2 ) 
3 (p,4>) 


ftr  • 


0 1 


at  these  points.  The  Implicit  Function  Theorem  implies  there  are 
c q > 0 and  a > 1 such  that  Equations  (2.8)  have  unique  solutions 
near  the  points  (p,4>)  above  for  |a|  < £Q,  |p|  + 1 6 1 < 1/a. 

Therefore,  in  the  original  variables  (p^,0,P4),  no  bifurcation 
occurs  in  the  region 


a,  e, 


= ( (P1»a,P4) 


S e0' 


a^/a} 


m 


This  implies  that  the  bifurcation  surfaces  in  an  neighborhood 

of  (p^fCr,p^)  = 0 must  be  in  the  region 

Ra,eQ  = {(-Pi'a'P4):  Pp  + P4  < Eo'  °2  - 01  ( I pi  I 3 + I P4  I 3 > > - 
In  the  region 


R1  = 


a,e, 


H { (P1»a,p4) : P4  > 0,  | p1 | < p41 


the  scaling 


2 2 3 

r = |Jpf  P4  = » Px  = f £7  = V|J  , n E 1-1,1],  )v|  < a(l  + J n| ) 

leads  to  the  new  bifurcation  equations 

— 3 3 2 

G,  = np  - j p + v cos  tj)  + pOCy  ) = 0 

2 (3.1) 

G2  = P + v sin  4>  + p0(y  ) = 0. 

For  y = 0,  one  obtains  directly  from  Equation  (3.1)  that 


a def 
A = 


det 


8 (GlfG2) 
8 (PA) 


1 


for  any  solution  (p,<|>)  of  (3.1). 


Since 


the  only  solution 


of  A = 0 is  p = 0.  But  this  implies  v = 0 which  implies  the 
Equation  (2.3)  is  autonomous.  For  this  situation,  we  know  Equation 
(2.3)  has  only  the  solution  v = 0 for  any  p.  > 0.  Therefore,  no 


bifurcation  occurs  in  R^. 

One  can  make  a similar  analysis  in  the  region 


R n 


{ (P1/Cf ,P4)  : P4  < 0,  |p1|  < -P4)  and  arrive  at  the  same 


conclusion . 


Consider  now  the  region 


R0  = R 
2 a,£. 


(P1fd,p4):  px  < 0,  | P4 1 < -P-,^}. 


2 3 4 _ 2 

If  we  use  the  scaling  r = y p , p-^  - -y  , cr  - vy  , p - my  , 
|m|  < 1,  | v | < a(l  + | m | ) the  bifurcation  equations  become 


— 3 3 2 

G1  = -p  - | p + v cos  <j>  + 0(y  ) 

_ 2 
G2  = mp  + v sin  <j>  + 0 (y  ) 


and,  at  y = 


3 (G^Gj) 
3 d , P ) 


p [- Cl  + I m2)  Cl  + I P2)  - rn2] 


This  quantity  is  zero  only  if  p = 0.  With  an  argument 
similar  to  the  one  for  R^  one  observes  there  is  no  bifurcation 
at  this  point. 

Let  us  now  analyze  the  region 


U = R H {(p1,a,p4):  p.  > 0,  | p4 | < PL). 

O U f £ q x ^ x 


For  the  scaling 


2 3 2 

r = yp,  Px  = y , a = vy  , p4  = my  , |m|  < 1,  |v|  < a(l  + |m|) 


the  bifurcation  equations  become 


— 3 3 2 

Gj^sp-^-p  + v cos  (j)  + pO(y  ) = 0 


= v sin  <p  + mp  + mpO(y  ) = 0. 


(3.2) 


Recall  that  we  need  only  consider  these  equations  for  |m|  < 1, 

| v | < a (1  + | m | ) and  p small.  Therefore,  p will  be  in  a 

1/3 

bounded  set,  |p|  < 2c  (1+a  ' ).  The  corresponding  determinant  of 

the  Jacobian  matrix  is  given  by 


A = det 


3 (G^ , G2 ) 

3 (P,<|>) 


q 2 2 

1 - \ P + 0(yZ) 


m + m0(y  ) 


- v sin  4) 


v cos  <p 


19  21  2 

= v 1 - j p cos  <J>  + m sin  <J>  + 0(y  ) 


For  v / 0 (i.e.  a / 0) , the  multiple  solutions  of  our  scaled 

bifurcation  equations  are  the  solutions  of  the  three  equations 


h-^  d=f  p - p3  + v cos  ()»  = p0(y2) 


J £ r\ 

h2  = mp  + v sin  <p  = mp0(£  ) 


def  Q j 2 

h^  = (1  - j p )cos  <J)  + m sin  <J>  = 0 (y  ). 


(3.3) 
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The  outline  of  the  program  is  now  as  follows. 

Treating  p as  a parameter,  we  determine  $,v,m  as  functions 
of  p such  that  Equations  (3.3)  for  y = 0 are  satisfied.  If 
det  d (hjyh^f h-j)/3  (<J>,v,m)  / 0,  then  we  can  determine  solutions  of 
Equations  (3.3)  for  y sufficiently  small.  Elimination  of  the 
parameter  p from  the  resulting  solutions  v*(p,y),  m*(p,y)  will 
give  the  possible  bifurcation  surfaces  v as  a function  of  m, 
or  perhaps  m as  a function  of  v,y.  We  then  verify  that  the 
number  of  solutions  (p,<J>)  of  Equations  (3.2)  [or  equivalently, 
the  number  of  solutions  (r,c)>)  of  Equation  (2.8)]  changes  by  two 
as  this  surface  is  crossed.  This  proves  that  it  is  a bifurcation 
surface . 

If  h^  = h2  = h^  = 0,  y = 0 then 


and  so 


also,  h.^  = h2  = 0 further  implies 


2 


v 


2 2 
m p 


(3.4) 


(3.5) 


(3.6) 


Equations  (3.4),  (3.5),  (3.6)  are  the  parametric  forms  of  the 
solutions  <(>,m,v  as  functions  of  p if  we  substitute  Equation  (3.5) 
into  Equations  (3.4)  and  (3.6).  Eliminating  p from  Equation  (3.5), 


(3.6),  we  have 


v2  = [i  + 9m2  ± (l-3m2)3/2].  (3.7) 

The  locus  of  the  points  in  the  (v,m)  plane  described  by  Equation 
(3.7)  are  plotted  in  Figure  1. 

If 

A-l  = det  9 (hlfh2,h3)/9  (<|>,m,v)  , 

then 

(9  2)  . 2 

A.  = v sin  <J>  - p x P - 1 cos  4>  sm  <p  - mp  cos  <J) . 

As  remarked  earlier,  if  A^  / 0 at  a solution  (<j>g,mg,  Vg,  pg)  of 
h^  = h2  = h^  = 0,  then  we  can  obtain  a solution  (4>,m,v)  of 
Equations  (3.3)  for  |p-pQ|  and  |y|  sufficiently  small  by 
employing  the  Implicit  Function  Theorem. 

If  h^  = h2  = h3  = 0,  then 

A^  = -2m 

2 

and  A^  = 0 implies  m = 0 since  m < 1.  If  m = 0,  then 
(v2,p2,4»)  = (16/81,  4/9,  0)  or  (0,  4/3,  n/2)  . 

Therefore,  at  each  solution  (<j>g,mg,  vQ,  pQ)  of  hl  = h2  = h3  = 
for  which  (Vg,mg)  / (±4/9,0)  or  (0,0),  there  is  an  e > 0 and 
unique  solutions  <f>  (p , y ) ,m (p,  y)  v (p , y ) of  Equations  (3.3)  for 
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|p“P0lf  |y|  < e,  <)>(p0,0)  = pQ,  m(p0,0)  = m0,  v(pQ,0)  = vQ.  Thus, 
if  we  exclude  a neighborhood  V of  the  points  (±4/9,0)  and  (0,0) 
in  the  (v,m)  plane,  then  we  can  find  a yQ  > 0 such  that  the 
Equations  (3.3)  can  be  uniquely  solved  for  (v,m,v)  as  functions 
of  (p,y)  for  |y|  < Pq  and  |p|  < 2c  (l+a'*'^ ) , the  a priori  bound 
on  p.  Elimination  of  p from  the  corresponding  functions  m,v 
gives  part  of  the  surface  indicated  in  Figure  2 in  terms  of  the 
unsealed  variables  (P1fCr,p4)  . The  first  approximation  to  the 
explicit  formula  for  this  surface  is  obtained  from  the  scaling  and 
Equation  (3.7)  and  is  given  by 


a = 


8 

81 


r ^ 

Lpi 


9PiP4 


± Ip?  - 


3P^>3/2 


2 


(3.8) 


At  the  points  (v,m)  = (±4/9,0)  in  the  (v,m)  plane  where 
A1  = 0,  we  compute  another  Jacobian  and  apply  the  implicit  function 
theorem.  At  these  points  and  for  (p,4>)  = (±2/3,0),  we  see  that 


A 


2 


det 


9 (h2 ,h2 ,h3) 

9 (p,m, v) 


= -2  / 0. 


Therefore,  we  can  solve  the  Equations  (3.3)  for  p,m,v  as  functions 
of  <f>,y  for  i{),y  sufficiently  small.  Eliminating  the  parameter 
<f>  gives  the  bifurcation  surface  near  the  points  (±4/9,0). 

At  the  point  (v,m)  = (0,0),  the  corresponding  solutions  (p,<)>) 
of  hx  = h2  = h3  = 0 is  (±2//3,tt/2)  . If  we  let  v = am,  then  the 
bifurcation  equations  are 


✓ 
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3 3 2 

= p - j p + am  cos  <f>  = 0(y  ) 

^2  . 2 
— = p + a sm  <p  = 0 (y  ) . 


(3.9) 


For  (p/<)>)  either  of  the  above  points,  we  have  h^/m  = 0 
implies  a = ±2//3.  We  wish  to  determine  the  multiple  solutions 
of  Equation  (3.9).  The  analysis  proceeds  as  before.  If 


A = det  3 (h1,h2/m)/3  (p,<)>) 


then 


~ 2 

A = a(l-9p  /4)cos  <J>  + am  sin  4> . 

We  now  consider  the  equations  h^  = 0,  h2/m  = 0,  A = 0 as  defining 
functions  p,m,a  as  functions  of  <p.  Along  the  solution  of  these 
equations,  we  have 

3 (h. ,h„/m,A)  . 

det  = ± ^ 0. 

3(p,m,a)  3/3 

Therefore,  we  can  find  functions  m((f>,y),  ot  ( , y ) and  p(<J>,y) 
that  will  satisfy  the  Equations  (3.9)  for  y close  to  zero  and  <!> 
close  to  ir/2.  Eliminating  $ from  the  functions  m(4>,y)  and 
a(<f>,y)  completes  the  discussion  of  the  possible  bifurcation  surface 
shown  in  Figure  2. 

It  remains  to  show  that  the  surface  in  Figure  2 is  the  bifurcation 


surface,  that  is,  we  must  show  that  the  number  of  solutions  ( p , <t> ) 
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I 


of  Equations  C2.8)  changes  as  this  surface  is  crossed.  However, 
we  need  not  check  all  points  on  the  surface.  It  is  sufficient, 
for  example,  to  show  that  the  number  of  solutions  changes  as  we 
cross  the  surface  in  the  plane  p4  = 0 since  the  only  possible 
way  for  the  number  of  solutions  to  change  is  to  pass  through  a 
multiple  solution.  For  p4  = 0,  Corollary  2.1  implies  all  solu- 
tions for  a = 0,  > 0 and  one  solution  for  p^  = 0,  o / 0. 

Therefore,  the  surface  in  Figure  2 is  a bifurcation  surface  with 
the  number  of  solutions  as  indicated. 

For  P2  < 0,  one  can  make  a similar  analysis  as  above  to  obtain 
the  corresponding  bifurcation  surface  lying  below  the  (cr,p4) 
plane.  Of  course,  one  has  two  more  periodic  solutions  in  each 
region  as  a result  of  Lemma  2.2.  We  do  not  draw  the  surface  in 
this  case. 

We  summarize  these  results  in  the  following  theorem. 

1/2  3 

Theorem  3.1.  Let  a = P2  p^.  There  is  a neighborhood  U in 

of  (0,0,0)  such  that  the  bifurcation  surface  T for  Equation  (1.1) 

with  (p1,a,p4)  E U is  depicted  in  Figure  2 and  the  surface  for 

P2  > 0 is  approximately  given  by  Equation  (3.8).  The  number  of 

2TT-periodic  solutions  of  Equation  (1.1)  at  a point  (plfa,p4)  c U 

which  are  continuous  in  (p3,p4)  is  shown  in  Figure  2. 

Remark  3.1.  The  results  of  this  paper  are  easily  extended  to  the 
equation 

77  + u - P]U  + P4  ar  - p2u3  + p3f  (tl 

at 


r 


provided  that  f (t)  is  an  even  function  of  t and 
•2tt 

f(t)cost  dt  = 7r.  One  uses  the  same  type  of  analysis  together 

0 

with  Remark  2.1. 
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